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ABSTRACT 

We have studied whether the jet in a collapse-driven supernova can be a key 
process for the rapid-process (r-process) nucleosynthesis. We have examined the 
features of a steady, subsonic, and rigidly rotating jet in which the centrifugal 
force is balanced by the magnetic force. As for the models in which the magnetic 
field is weak and angular velocity is small, we found that the r-process does not 
occur because the final temperature is kept to be too high and the dynamical 
timescale becomes too long when the neutrino luminosities are set to be high. 
Even if the luminosities of the neutrinos are set to be low, which results in the 
low final temperature, we found that the models do not give a required condition 
to produce the r-process matter. Furthermore, the amount of the mass outflow 
seems to be too little to explain the solar-system abundance ratio in such 
low-luminosity models. As for the models in which the magnetic field is strong 
and angular velocity is large, we found that the entropy per baryon becomes 
too small and the dynamical timescale becomes too long. This tendency is, of 
course, a bad one for the production of the r-process nuclei. As a conclusion, we 
have to say that it is difficult to cause a successful r-process nucleosynthesis in 
the jet models in this study. 

Subject headings: nucleosynthesis, abundances — stars: magnetic — stars: 
rotation — supernovae: general 
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1. INTRODUCTION 



It is one of the most important astrophysical problems that the sites where the 
rapid-process (r-process) nucleosyntesis occurs are not still known exactly. There are, at 
least, three reasons that make the study on r-process nucleosynthesis important. One of 
them is a very pure scientific interest. The mass numbers of the products of r-process 
nucleosynthesis are very high (A = 80-250), which means that the most massive nuclei in 
the universe are synthesized through the r-process. You can guess easily that the situation 
in which the r-process nuclei are synthesized is a very peculiar one in the universe. We want 
to know where, when, and how the r-process nuclei are formed. Second reason is that some 
r-process nuclei can be used as chronometers. For example, the half-lives of 232 Th and 238 U 
are 1.405 xlO 10 yr and 4.468 xlO 9 yr, respectively. So if we can predict the mass-spectrum of 
the products of r-process nucleosynthesis precisely, we can estimate the ages of metal-poor 
objects which contain the r-process nuclei by observing its abundance ratio. Third reason 
is that some r-process nuclei can be used as tools of the study on the chemical evolution in 
our Galaxy (e.g., |Ishimaru and Wanajo 1999|) , which has a potential to reveal the history 



of the evolution of our Galaxy itself. Due to the reasons mentioned above, the study on the 
r-process nucleosynthesis is very important. 

The conditions in which the r-process nucleosynthesis occurs successfully are (e.g., 
Hoffman et al. 1997] ): (i) neutron-rich (n n > 10 20 cm -3 ), (ii) high entropy per baryon, 



(hi) small dynamical timescale, and (iv) small Y e . This is because r-process nuclei are 
synthesized throuth the non-equilibrium process of the rapid neutron capture on the seed 
nuclei (iron-group elements). In other words, a highly dense, explosive, and neutron-rich site 
with high entropy will be a candidate for the location where the r-process nucleosynthesis 
occurs. 
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The most probable candidates for the sites are collapse-driven supernovae (e.g., |Wooslcy 



et al. 1994) and/or neutron star mergers (e.g., |Freiburghaus et al. 1999]) . This is because 



these candidates are thought to have a potential to meet the requirements mentioned above. 
However, we think that the collapse-driven supernovae are thought to be more probable 
sites than the neutron star mergers, because metal poor stars already contain the r-process 
nuclei (e.g., I-Kreiburghaus et al. 1999| ). In fact, Mcwilliam et al. (1995) reported that the 



abundance of Eu can be estimated in 11 stars out of 33 metal-poor stars. These observations 
prove that r-process nuclei are produced from the early stage of the star formation in our 
Galaxy. Taking the event rate of collapse-driven supernovae (10~ 2 yr _1 Gal -1 ; [van den 



Pergh and Tammann 199~T|) and neutron star merger (10 5 yr 1 Gal 1 ; |van den Heuvel and 



Lorimer 1996| ; Pehte and Brown 1998|) into consideration, collapse-driven supernovae are 



favored since they can supply the r-process nuclei from the early stage of the star formation 
in our Galaxy. Also, Cowan et al (1999) reported that the abundance ratio of r-process 
nuclei in metal poor stars are very similar to that in the solar system. This proves that 
r-process nuclei are synthesized through the similar conditions. This will be translated that, 
at least, most of the r-process nuclei are from one candidate. So we assume in this paper 
that most of the r-process nuclei are synthesized in the collapse-driven supernovae. 

There are many excellent and precise analytic and/or numerical calculations on the 
r-process nucleosynthesis in the collapse-driven supernovae. However, it would be able to 
be said that there is no report that the r-process nuclei can be reproduced completely. For 
example, Takahashi et al. (1994) performed numerical simulations assuming Newtonian 
gravity and reported that entropy per baryon in the hot bubble is about 5 times smaller 
than the required value. Qian and Woosley (1996; hereafter QW96) also reported analytic 
treatments of the neutrino-driven winds from the surface of the proto-neutron star. At the 
same time, their analytical treatments are tested and confirmed by numerical methods. 
However, the derived entropy by their wind solutions are shown to fall short, by a factor 
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of 2-3, of the value required to produce a strong r-process ( [Hoffman et al. 1997|) . Otsuki 
et al. (2000) have surveyed the effects of general relativity parametrically. They reported 
that r-process can occur in the strong neutrino-driven winds {L v ~ 10 52 erg s _1 ) as long 
as a massive (~2.0 M & ) and compact (~ 10 km) proto-neutron star is formed. It is very 
interesting because such a solution can not be found in the frame work of Newtonian 
gravity flQian and Woosley 1996|). Such a solution is confirmed by the excellent numerical 



calculations ( |Sumiyoshi et al. 1999j ). However, the equation of state (EOS) of the nuclear 



matter has to be too soft to achieve such conditions. In fact, the r-process nuclei can not 
be produced in the numerical simulations with an normal EOS ( [Sumiyoshi et al. 1999Q . 

There is only one report that r-process nucleosynthesis occurred successfully. That 
is the work done by Woosley et al (1994; here after WWMHM94). In their numerical 
simulation, the entropy per baryon becomes higher and higher as the time goes on. Finally, 
at very late phase of neutrino-driven wind (~ 10 s after the core-collapse), successful 
r-process occurs. However, there are some questions on their results. One is that the reason 
why the entropy per baryon becomes high at the late phase is unclear. In fact, when we 
adopt the analytic formulation of QW96, it is concluded that such a high entropy should 
not be obtained. It is true that general relativistic effects are included in WWMHM94, 
but such a high entropy could not be obtained in the work done by Otsuki et al (2000). 
So the discrepancy between WWMHM94 and QW96 can not be explained by the general 
relativistic effects alone. Also, WWMHM94 has a problem that the nuclei whose mass 
numbers are ~ 90 are much produced in the early stage of the neutrino-driven winds. 
So if we try to reproduce the mass-spectrum of the solar system abundances, we have 
to abandon the matter which is made at the early stage of the neutrino-driven winds. 
Even worse, the successful mass-spectrum at the late phase of neutrino-driven winds is 
destroyed when the reactions of neutral-current neutrino spallations of nucleons from 4 He 
are taken into consideration QMeyer 1995|) . In their paper, it is reported that (30-50)% 
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increase in the entropy is needed in order to restore the A = 195 peak, which is extremely 
large modifications to the model. So, although the work done by WWMHM94 is the very 
remarkable and interesting one, it would not be concluded that the problem of r-process 
nucleosynthesis has been solved completely. 

Due to the reason mentioned above, it will be natural to think that there may be 
an (some) effect (s) that will help the r-process nucleosynthesis. There are some papers 
that mention the effects of jet ( (Symbalisty 1984] ; |Shimizu, Yamada, k Sato 1994j ; [Nagatakj| 
[2000|) that is generated around the polar region of the highly rotating (the period T ~lms) 
and/or magnetized (~ 10 15 G) proto-neutron star ( [LeBlanc k Wilson 197U| ; |Symbalisty| 



|Shimizu, Yamada, k Sato 1994] |Yimada k Sato HgM|; |Fryer k Heger l999| ). Since the 



physical conditions in the jet are different from those in a spherical explosion, the products 
of the nucleosynthesis are expected to be different. In fact, it is reported that the products 
of the explosive nucleosynthesis in a jet-like explosion in Si- and O rich layer are much 
different from those in a spherical explosion ( [Nagataki et al. 1997| ; |Nagataki 2000| ). So it 



is natural to think that the products of explosive nucleosynthesis in the Fe core (that is, in 
the hot bubble) will be also changed due to the effects of the jet. This is the reason why we 
should examine whether the effects of the jet can help the synthesis of r-process nuclei or 
not. 

Since there is no numerical simulation on the r-process nucleosynthesis in the jet 
during the neutrino-driven wind phase, our final goal is to perform such realistic numerical 
simulations. However, such a numerical simulation will be a heavy task. We also think 
that the results of such numerical simulations will not be understood completely without 
analytical studies. So, before performing such numerical simulations, we examine the 
physical conditions of the jet using a simple model. In this paper, we study a steady, 
subsonic, and rigidly rotating flow of the jet that is driven by neutrinos. We study the 
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effects of the jet on the physical conditions such as the entropy per baryon and the 
dynamical timescale. Finally, we discuss whether the effects of the jet can help the synthesis 
of r-process nuclei. 

In section |2], we explain the formulation for the jet in the hot bubble. In section |3|, we 
show the results. Summary and discussions are presented in section [|. 

2. FORMULATION OF JET 
2.1. Basic Equations 



In Gaussian units, the Euler's equation acted on by electromagnetic forces can be 
written as ( [Shapiro and Teukolsky 1S5B| ) 

v ^ = --VP- V$-— VB 2 + — (B- V)B. (1) 

at p 8irp Aixp 

Here 

d 9 _ _ . . 

is the Lagrangian time derivative following a fluid element. 

In this paper, we study a steady jet which has 0-symmetry around the polar region of 
the proto-neutron star. So we use the cylindrical coordinate (r, 0, z) for convenience. The 
origin z=0 is set at the center of the proto-neutron star. In the cylindrical coordinate, the 
Euler equation can be written as 

/ ftv , „ ftv £ \ / 8P x / n v / (OBz-BBAb,- ^HrB^) \ 



O \ dP \ QM 

X dz ' ' - ' 



(3) 
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Conservation of mass requires that 



rD 

M = 2irrpv z dr, 
Jo 



(4) 



where M and D are the constant mass outflow rate in the ejecta and the radius of the jet, 
respectively. 

The equation for the evolution of material energy, e, is 



PQ 



V ■ (pev) + PV ■ v 

f 1 9 / \ d . \ /IS. . 9 



(5) 
(6) 



where q is the net specific heating rate due to neutrino interactions ( |Qian and WoosIey| 
|1996| ). In this study, we consider three neutrino heating and/or cooling processes (neutrino 
absorption on free nucleons, neutrino scattering processes on the electrons and positrons, 
and electron capture on free nucleons) as 



Q — QuN + Que — QeNi 



(7) 



where 



q vN = 1.55 x 10' 5 N A 



1 Y e )L Ve ^ie 2 u MeV + ^e-^fe,51 e L,MeV 



1 — X 

Rr 



erg s 1 g 1 



q ve = 3.48 x KT'Wa 
and 



6 AT MeV i £ 



6 \1-X 

v e ,51 e ^ e ,MeV + -^P e ,51 e i/ e ,MeV + ^-^^,51 e ^,MeV ) — 

( J Kq 



erg s 1 g 1 



(9) 



q eN = 3.63 x 10- 6 N A T^ cY 



erg s 1 g 1 



(10) 



Here Rq is the neutrino sphere radius in units of 10 6 cm, x = (1 — R 2 /z 2 ) 1 ^ 2 , is 
Avogadro's number, L^ 5 i is the individual neutrino luminosity in 10 51 erg s _1 , e^MeV is an 
appropriate neutrino energy e u in MeV (|Qian and Woosley 1996| ). In this study, we set q=0 
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at T <0.5 MeV, because free nucleons are bound into a-particles and heavier nuclei and 
electron-positron pairs annihilate into photons. 

The pressure P and internal energy e are determined approximately by the relativistic 
electrons and positrons and photon radiation as long as T > 0.5 MeV. So the pressure and 
internal energy can be written as 



llvr 2 fc 4 T 4 
180 k 3 c 3 



and 



llvr 2 fc 4 T 4 



dyn cm 



erg g 



(11) 



(12) 



60 h 3 c 3 p 

where k and h are Boltzmann and Planck constants, respectively. These are the basic 
equations in this paper. Precisely, we should consider the effects of annihilations of 
electron-positron pairs on the dynamics at T < 0.5 MeV. However, its effects seems to be 
little flSumiyoshi et al. 1999| ). 



2.2. Model for the Jet 



In this study, we try to construct a simple steady state for a jet so that we can 
understand the effects of magnetic fields clearly. So we assume in this study that the form 
of the magnetic fields is 

3= (0,0, B+faz)). (13) 

Next, we seek a solution for the jet rotating rigidly. That is, we seek a solution v = 
(0,Qr,v z (r, z)) at z > R, where Q is the angular velocity of the proto-neutron star. From 
Eq. (|3p, B and v must meet the following relation: 



-10- 



In this case, there is a useful approximate solution. That is, 

fl* = (27rp)tor, (15) 

where p and P are assumed to be a function of z alone. It is noted that this approximate 
solution is valid as long as QD is small enough. In fact, the derivatives of p and T by z can 
be written as 



dT _ i + (f + $) 



(16) 



dp 

*~ $ (17) 

It is clearly understood from Eq. ( |I6|) that the solution holds valid as long as QD is small 
enough. In this case, Eq. (f|) can be rewritten as 

M = nD 2 pv z . (18) 



From Eqs. (^4|) and (|T^), we can see that the magnetic fields accelerate the matter as 
long as the density gradient is negative. We are also able to see clearly from Eq. (|i"5|) that 
small QD corresponds to weak magnetic field. So we call such solutions weak field solutions 
in this paper. On the other hand, we also want to see what can happen as the magnetic 
fields are getting stronger and stronger. So, although the truth is that these solutions are 
inapplicable to the case of strong magnetic fields, we extend these solutions to the case in 
which QD is relatively large. By so doing, we will be able to find at least what happens in 
the case of strong magnetic fields qualitatively. We call such solutions strong field solutions 
in this paper. In section we show the results for the weak and strong field solutions, 
respectively. 
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2.3. Boundary Conditions 



In this study, the surface of the proto-neutron star is considered as the inner boundary. 
The inner boundary conditions are composed of density, luminosities of neutrinos, mass and 
radius of the proto-neutron star, velocity of the outflow, and QD. Initial temperature and 
final electron fraction are determined by these parameters as ( |Qian and Woosley 1996|) 



Tj = 1.19 x 10 



10 



i i 



Ll M ^4 eM eV [K] (19) 



and 



Y = ( 1 + Lpe e ^,Mev - 2 A + 1.2 A 2 / e p<;i Mev \ 
\ L Ue ej, ej MeV + 2A + 1.2A 2 / e^MeV / 

where L v ^\ is the individual neutrino luminosity in 10 51 ergs s _1 , Rq is the neutron star 
radius in 10 6 cm, A = 1.293 MeV is the neutron-proton mass difference, and e^MeV is a 
neutrino energy in MeV. We assume that the neutron star radius is equal to the neutrino 
sphere radius. 



In this study, the luminosities of neutrinos are assumed to be common QQian and 



Woosley 1996| ; |Utsuki et al. 200"UD . The energy of neutrinos are assumed to be 12, 22, 



and 34 MeV for z/ e , u e , and other neutrinos, respectively ( |Woosley et al. 1994j ; |Cjian and 



Woosley 1996| ; ptsuki et al. 2000] ). Surface density is assumed to be 10 10 g cm 3 ( |Otsuki 



et al. 2000| ). Initial velocity of the outflow is chosen so that M becomes less than M, 



crit; 



where M crit is the critical value for supersonic solution. We have to emphasize that the flow 
is assumed to be subsonic and contain no critical point. Previous works also adopted this 
assumption QQian and Woosley 1996| ; Ptsuki et al. 2000| ). In case we try to survey the 



flow that contains a transition point like a shock front, we can not use Eqs. @, (|i~4l) , (fil|), 
and (|i~7|). This is because these differential equations diverge and break down. In order 
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to treat such a flow that contains a discontinuity, we have to use the Rankine-Hugoniot 
relation instead of these equations. We will examine such flows in the forthcoming paper. 
D is assumed to be 10 5 cm, which is about one tenth of the proto-neutron star radius. We 
note that an appropriate environment can exist around the jet, that is, there is a solution 
that meets Eq. (|I4 ) for r > D in which and becomes zero at r — > oo. So the rigid 



rotating flow at r < D is not an unrealistic solution at all. In this study, we consider only 
the region r < D in which the matter rotates rigidly for simplicity. Other parameters are 
changed parametrically. The parameters employed as well as the model names are given in 
Table [T| We take z = 10 9 cm for the radius of the outer boundary ( ptsuki et al. 20"00| ). 
Since the radius of the Fe core is about 10 8 cm, we think the radius of the outer boundary 
is large enough to investigate the r-process nucleosynthesis that occurs in the hot bubble. 

EDITOR: PLACE TABLE | HERE. 



3. RESULTS 



3.1. Weak Field Solutions 



Output parameters are shown in Table [| Entropy per baryon (S) , dynamical timescale 
(Tdyn), analytically estimated dynamical timescale (rd yn . a na), electron fraction (Y e ), and 
temperature (T) at the outer boundary (z = 10 9 cm) are shown in the table. Entropy per 
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baryon in radiation dominated gas can be written as 

b/ 45 sVp/m» ( ' 

as long as T > 0.5 MeV. is the nucleon rest mass. Even if the annihilation of the 

electron-positron pairs occurs, the entropy per baryon is conserved. The definition of the 

dynamical timescale is the time for the temperature to decrease from 0.5 MeV to 0.2 MeV. 

This is because r-process nucleosynthesis occurs in this temperature range ( |Woosley et al.| 

1994 ; lakahashi et al. 1994 ; |Qian and Woosley 1996| ). The reason why the dynamical 



timescales are not written in some models is that temperature does not decrease to 0.2 
MeV within z = 10 9 cm. The definition of the analytically estimated dynamical timescale 
is described later (see Eq. (|3ll) ). 

EDITOR: PLACE TABLE | HERE. 

There is a tendency that S and Tf, in the models studied here are higher than those in 
the models of QW96. For example, we show in Figure [I] the outflow velocity, temperature, 
and density as a function of z for Model 10CW and Model 10C in QW96. The inner 
boundary conditions in Model 10C in QW96 is same as those in Model 10CW (see table 
1 in |Qian and Woosley 1996 ). From this figure and Eq. ( pl[) , we can see the tendency 



mentioned above clearly. We consider the reason for this tendency below. In order to find 
a clue to it, we show the absolute values (in cgs units) of the components in Eqs. ( ^6|) 
and ( jl7|) for Model 10CW in Figure |2j. From this figure, we can see what determines the 
gradients of temperature and density. For comparison, we also show those for Model 10C in 
QW96 in Figure |. Lines (a)-(g) correspond to 4e/T, (P/ep)(e + P/p), v 2 z , pGM/z 2 , Pq/v z e, 
(e/p + P/p 2 )^, and q/v z as a function of z, respectively. As can be seen from Figures [2] 
and [| dp/dz can be approximated by 

dp 



dz 



V p) 



(22) 
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Also, dT/dz can be approximated by 



dT + m N GM 



dz ~j=r S z^ 



(23) 



From Eq. (|23| ) we can see that the temperature gradient dT/dz is steeper when the entropy 
per baryon is smaller. So we can expect that the entropy per baryon is larger in Model 
10CW than Model IOC in QW96. It is verified in Figure |. So, the problem why T b 
becomes high is translated to the problem why the entropy per baryon becomes high in 
Model 10CW. The answer is shown in Figure |2|. The velocity at small radii in Model 10CW 
is smaller than that in Model IOC in QW96. The entropy per baryon is getting larger with 
time as long as the heating process is effective. The heating process is effective only at 
small radii, which is proved in Figure f|. So, if the velocity at small radii is small, there is a 
long time for the outflow matter to get heat from neutrinos. As a result, higher entropy per 
baryon can be achieved and T& becomes higher. This is the case with the jet models studied 
in this paper. As for the reason why the initial velocity has to be smaller in Model 10CW 
than Model IOC in QW96 is relatively difficult. At least, we can say that dv z /dz can be 
written from Eq. fll8[) as 

<t L = _v 1 _dp inModdl0CW> (24) 
dz p dz 

On the other hand, it can be written as 

du 1 = _v ± d^_v 1 in Model 10C in QW96) (25) 
dz p dz r 



because 4:7rz 2 pv z is constant in Model IOC in QW96. Taking these equations and Eq. (|16|) , 
which also holds in the models in QW96 ( IQian and Woosley 1996| ), into consideration, the 
velocity gradient is steeper in Model 10CW as long as v z (z), p(z), and T(z) are same in 
both models. So we can guess that initial v z has to be smaller in Model 10CW so as not to 
achieve the adiabatic sound speed, v s = (4P/3p)^. This will be the reason why the initial 
velocity in Model 10CW has to be smaller. 
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EDITOR: PLACE FIGURE |T| HERE. 



EDITOR: PLACE FIGURE g HERE. 



EDITOR: PLACE FIGURE |3| HERE. 



EDITOR: PLACE FIGURE | HERE. 



Next, we consider the dynamical timescale. As mentioned above, the dynamical 
timescale becomes infinity when TJ, is above 0.2 MeV. Since TJ, in the jet models tends to 
be higher, we can not determine the dynamical timescale in many jet models. This is a 
bad tendency to produce r-process nuclei. From Table [§ the dynamical timescale can be 
determined as long as £p ei 5i is small. We will discuss later whether r-process can occur or 
not in such models. Before we continue to further discussion, we explain the way to derive 
the analytically estimated dynamical timescale rd yniana . 

Like Model 10C in QW96, from Eqs. (§) and (|14]) we can derive the following equation 
as long as QD is small: 



So, when the heating/cooling process becomes non-effective, there is a conserved quantum 




(26) 




(27) 



In the subcritical case, v z < v s everywhere, we can approximate Eq. fl2"7j ) as 



TS GM T b S 



(28) 



m N 



r 



m N 
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When T h -C0.5MeV, we can take 



TS GM 



m N 



at 0.5MeV. 



(29) 



See Qian and Woosley (1996) for details. In this case, T oc r _1 near T ~ 0.5 MeV because 
S is also a conserved quantum as long as the heating/cooling process does not work. From 
Eq. (PI), p oc r -3 . From Eq. fllSl) , v z oc r 3 . So, r dynana can be approximated as 



yn,ana 



T=0.2MeV /-T=0.2MeV dt dz 

dt = ——dT 

T=0.5MeV JT=0.5MeV dz dT 

r \ /-T=0.5MeV / r - 

— ) / TrfT- 0.42 

z / T=0.5MeV JT=0.2MeV \V / T=0.5MeV 



(30) 

M • (31) 



sec 



For further discussion, we try to express the dynamical timescale by input parameters 
Qian and Woosley 1996 for details). From simple calculations, we can derive the 



entropy per baryon and the mass outflow rate in Model 10CW as 

M \ 



S/k ~ 235L^ 51 e Pe \ leV _R 6 



I AM, 



0, 



and 



M 



D 2 
4R2 



x 1.14 x 10" 10 L 3 



in 

3 



, o D3 

P e ,51 e i/e,MeV- rt 6 



I (\AMr. 







M 



M & s- 1 



respectively. At last, we can derive the expression for rd yn ,ana from Eqs. 
(m), and (131) that 



r dyn,ana ~ ^-^Lpf^ep^^yRQ 



M 



[a]. 



(32) 



(33) 



(34) 



v 1.4M 0/ 

This value for each model is written in Table 0. We can find that Td yn ,ana agrees with 
Td yI1 within a factor of 2-3 in all models in which the dynamical timescale can be defined, 
although Td yI1 ,ana tends to be smaller. This precision of Td yni ana is about the same that in the 
analysis in QW96 (see |Qian and Woosley 1996| for details). 
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We compare r dyn ana derived here with that in QW96 in order to see the effects of jet. 
It is written as 

2 

Ilr^L = 0.15—^-^1 . (35) 

r dyn , Q w Ll ej51 eI ejMeV 

For example, when Lp ej5 i = 0.1, R 6 = 1, and e PeiM eV = 22, r dyn /r dyniQW = 0.04. Even if 
we take the uncertainty (a factor of 2-3) of 7d ynjana into consideration, it seems that short 
dynamical timescale can be realized in jet models. So, as long as T& can be small (that 
is, Lj? e) 5i is small enough), we think that short dynamical timescale can be realized in jet 
models. 

Now we discuss whether r-process nucleosynthesis can occur in jet models or not. We 
can use the Hoffman's criterion ( [Hoffman et al. 1997| ) for the judgment. The Hoffman's 
criterion can be written as 

i 

S > 2 x 10 3 r e f^V (36) 



for Y e >0.38. This is the criterion for production of r-process nuclei with mass number 
A ~200 QHoffman et al. 1997Q . Substituting Eqs. (0) and fl3§ into Eq. We can 



translate the criterion into 

0.74 a / M \" 
i? 6 <^L? 2 51 — — , (37) 

where ep ei M e v is assumed to be 22. C is the factor 2-3 which represents the tendency that 
Tdyn.ana tends to be smaller than Td yn . In case of M=1.4M and C=2, we can see from 
Eq. ( PT| ) that Rq can be larger than 1 as long as L Pe 5 i is larger than 8.6. This means that 
there seems to be possible to produce r-process nuclei without unrealistically soft EOS 
as long as Lp e ,5i is larger than 8.6. However, we found from Table |2] that Tb can not be 
lower than 0.2 MeV when the neutrino luminosity is so high. So we have to conclude that 
r-process can not occur in these models. 
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However, we can find that the requirement to produce r-process nuclei is relaxed in 
these jet models. In fact, for the models of QW96, the criterion like Eq. (|37|) can be written 

as 

R,,qw < 0.19L5 ei51 [jj^J ■ (38) 

So, i?6,QW can be larger than 1 when L Pej51 is larger than 2.3 xlO 4 , which can not be realized 
in numerical simulations of supernovae flWoosley et al. 1994|) . We can compare the required 
radius in jet models (C = 2) with that in models of QW96 as 

Rr 2 ( M \ 33 

-^- = 3.2Lf 51 . (39) 

For example, in case of M — 1.4M Q and L Pe 5 i = 0.1, -R 6 /-R 6) qw = 2.8. This suggests that 
the EOS of the nuclear matter has not necessary to be too soft as Otsuki et al (2000) 
requires, when the jet models are adopted. So, we can conclude that the requirement 
to achieve r-process nucleosynthesis is relaxed for the jet models, although r-process 
nucleosynthesis is not still able to occur in the jet models. When we include the effects of a 
jet and general relativity at the same time, an successful solution may be obtained. We will 
discuss this topic in the forthcoming paper. 

Finally, we consider the mass outflow rate. Taking the event rate of collapse- driven 
supernovae (10~ 2 yr _1 Gal -1 ; |van den Bergh and Tammann 1991| ), one needs 10~ 6 to 



10 4 Mq of ejected r-process material per collapse-driven supernova event to explain the 
observed solar-abundance ratio QKappeler et al. 1989| ; |Woosley et al. 1994|) . We give a 



rough estimate whether each model can produce r-process matter enough to explain the 
solar-system abundance. Taking the gravitational binding energy of a neutron star is 
~ 3 x 10 53 erg into consideration, the duration time, t, can be roughly estimated as 
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where the factor 6 represents the number of flavor of neutrinos. So required M req to explain 
the solar abundance ratio is 

M req > (l(T 4 - 1(T 6 ) M x ^1 [m & s' 1 } . (41) 

So M in Table |I| has to meet the relation 

2xi~ M req (42) 

where the factor 2 represents that the matter is ejected from both (north and south) sides of 
the polar region. We can find from Eq. (f|l]), (0), and Table [l] that low- luminosity models 
like 10GW-10LW and 30FW-30HW can not meet the relation mentioned above, although 
the high-luminosity models like 10AW-10FW and 30AW-30EW may be able to meet the 
relation. 

At last we can derive the conclusion. In the framework of the steady neutrino-driven jet 
with weak magnetic fields, high-luminosity models can not work because T& becomes high 
and dynamical timescale becomes too long. Also, low-luminosity models do not present a 
required condition to produce the r-process nuclei in the framework of Newtonian gravity. 
Although a successful solution may be obtained when we include the effects of jet and 
general relativity at the same time, the amount of mass outflow seems to be too little to 
explain the solar-system abundance ratio in such low-luminosity models. As a conclusion, 
we have to say that it is difficult to cause a successful r-process nucleosynthesis in the weak 
field solutions. 



3.2. Strong Field Solutions 



In this subsection, we examine the models in which QD ^0. As stated in subsection 771 
the truth is that the solutions studied in this paper are inapplicable to the case of strong 
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magnetic fields. However, we want to see what can happen as the magnetic fields are getting 
stronger and stronger. At the same time, we should investigate the range of application of 
weak field solutions. 

In this study, we consider three models (Model 10GS1-3) for strong field solutions. 
These are listed in Table [I]. The strength of the magnetic fields corresponds to 2.5 xlO 13 G, 
2.5xl0 14 G, and 2.5xl0 15 G, respectively (see Eq. (P)). 

The results are shown in Table [| We can find that the values of output parameters 
change drastically in Model 10GS3. This result suggests that the weak field solutions can be 
applied at least in the range of order B < 10 14 G. This means that the weak field solutions 
will be valid for the typical proto-neutron stars. In case we consider the nucleosynthesis in 
a magnetar ( Woods et al. 1999| ), we will have to investigate using strong field solutions. 



In order to understand what happens in Model 10GS3, we show the outflow velocity, 
temperature, density, and entropy per baryon as a function of z of Model 10GS3 in Figure |[ 
We also show the absolute values (in cgs units) of the components of Eqs. ([R]) and ( |T7D for 
Model 10GS3 in Figure §. Lines (a)-(h) correspond to 4e/T, (P/ep)(e + P/p), v 2 , pGM/z 2 , 
Pq/v z e, (e/p + P/p 2 )^, q/v z , and VL 2 D 2 /4 as a function of z, respectively. As can be seen 
from Figures ||, dp/dz at r = D can be approximated by 



dp 

dz 



pGM/z 2 



(43) 



Also, dT/dz can be approximated by 

dT T(z + %)% 
dz 4e 

Tq 



4ev, 



for z < 8 x 10 1 and z> 2 x 10 3 (44) 
for 8 x 10 1 < z < 2 x 10 3 . (45) 



EDITOR: PLACE FIGURE | HERE. 
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EDITOR: PLACE FIGURE | HERE. 

We note that the denominator of Eq. (|43| ) is Q 2 D 2 /4, which is constant through out. 
Since the denominator is large, the value of the density gradient becomes small. So the 
density gradient in 10GS3 becomes nearly zero at smaller z (compare Figures |T] and [JJ). 
As a result, temperature does not decrease due to the adiabatic expansion. When the 
temperature is kept to be high, cooling process works better (see Eqs. (^), (Q), and (|I0D). 
So, at z ~ 8 x 10 6 cm, the cooling process overcomes the heating process. As a result, 
entropy per baryon is getting lower in the range 8 x 10 6 < z < 2 x 10 8 cm. This is the 
reason why the entropy per baryon becomes small in Model 10GS3. We also note that the 
temperature gradient becomes too small at z > 2 x 10 s cm. This is because the numerator 
of Eq. (44) becomes in proportion to dp/dz, which is very small as seen above. So we can 
find that the entropy per baryon is small and the dynamical timescale is long in the strong 
field solution 10GS3. This tendency is a bad one for the production of r-process nuclei. The 
truth is, of course, that the solutions studied in this paper are inapplicable to the case of 
strong magnetic fields (see Eq. (|j~6D). However, similar tendency discussed here may be also 
found in the exact solutions. If so, the strong field solutions are unsuitable for the r-process 
nucleosynthesis. 

Finally we add a comment on the effects of the magnetic fields. You may have thought 
that the dynamical timescale in the strong field solutions will become shorter than that in 
the weak field solutions, because the magnetic fields accelerate the matter as long as the 
density gradient is negative (see Eqs. (|T4]) and (|l~5|)). However, in the framework of the 
steady flow, the density gradient becomes nearly zero and acceleration does not occur. So, 
the situation may be different and acceleration occurs in the case of an unsteady flow. We 
will investigate the features of such an unsteady flow by numerical tests in the near future. 
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4. SUMMARY AND DISCUSSION 



We have studied whether the jet in a collapse-driven supernova can be a key process 
for the r-process nucleosynthesis. We have studied the effects of a jet using a simple model, 
because the results of a realistic numerical simulation concerning with such a theme will not 
be understood completely without analytical study. Although our final goal is to perform 
such realistic numerical simulations, this work is a necessary process to understand the 
effects of a jet on the r-process nucleosynthesis. 

We have studied two cases, that is, the cases in which the magnetic fields are 
weak/strong. In both cases, we assumed that the flow is steady and subsonic. 

As for the weak field solutions, we have concluded that high-luminosity models can 
not work because T b becomes high and dynamical timescale becomes too long. Also, 
low-luminosity models do not present a required condition to produce the r-process nuclei 
in the framework of Newtonian gravity. Although a successful solution may be obtained 
when we include the effects of jet and general relativity at the same time, the total mass 
of the outflow seems to be too little to explain the solar-system abundance ratio in such 
low-luminosity models. On the other hand, we found that the entropy per baryon is small 
and the dynamical timescale is long in the strong field solution. We can tell that this 
tendency is a bad one for the production of r-process nuclei. The truth is that the solutions 
studied in this paper are inapplicable to the case of strong magnetic fields. However, similar 
tendency may be also found in the exact solutions. If so, we will be able to conclude that 
strong field solutions are unsuitable for r-process nucleosynthesis. As a conclusion, we have 
to say that it seems to be difficult to cause a successful r-process nucleosynthesis in the jet 
models in this study. 
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We have to emphasize that there are some assumptions in this study. So, we can not 
say that we have proved that a successful r-process nucleosynthesis does not occur in a 
neutrino-driven jet in a collapse-driven supernova. For example, we assumed that the flow 
is subsonic and there is no critical point, which is the common assumption in the previous 
studies on the r-process nucleosynthesis ( |Qian and Woosley 1996| ; |Otsuki et al. 2000|) . 
However, we think that we do not need to restrict the solutions in such a way, that is, there 
may be a transition point at which Eqs. (Q), (13), ([H]), and ( |T7j ) break down. In most cases, 
the transition point will be a shock front. It means that the flow will gain entropy at the 
transition point, which will be a good sense to produce the r-process nuclei. The problem 
whether the flow contains transition points or not depends sensitively on the initial velocity 
on the surface of the proto-neutron star. So our final goal is to determine physically the 
velocity at the surface of the proto-neutron star. It means that the M should be not given 
as an input parameter. It should be an output parameter. We have to investigate the 
mechanism for determining the outflow velocity at the surface of the proto-neutron star for 
further discussions. We also assumed that the flow is steady. As seen in subsection CO, the 
density gradient becomes nearly zero and acceleration due to the strong magnetic fields 
does not occur in the framework of the steady flow. So, the situation may be different and 
acceleration occurs in the case of the unsteady flow. We should investigate the features 
of the unsteady jet flow. It will be investigated by numerical tests assuming a simple 
environment. We will perform such numerical tests in the near future. At the same time, 
we will try to perform realistic numerical simulations for the jet induced explosion in a 
collapse-driven supernova and seek a final answer whether the jet in a collapse-driven 
supernova is a key process for the r-process nucleosynthesis or not. 
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Fig. 1. — Outflow velocity, temperaure, and density as a function of height (z) from the 
center of the proto- neutron star. These quanta are written in unit of 10 7 cm s _1 , 1 MeV, 
and 10 8 g cm -3 , respectively. Solid lines correspond to Model 10CW, whereas dashed lines 
correspond to Model IOC in QW96. 
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Fig. 2.— Absolute values (in cgs units) of the components of Eqs. (16) and (17) for 
Model 10CW. Lines (a)-(g) correspond to 4e/T, (P/ep)(e + P/p), v 2 z , pGM/z 2 , Pq/v z e, 
(e/ p + Pj p 2 )^f, and q/v z as a function of z, respectively. The discontinuities of lines (e) and 
(g) at z ~ 10 7 cm reflect the freezeout of the neutrino reactions (see subsection 2.1). 
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Fig. 3.— Same as Figure 2, but for Model IOC in QW96. 
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Fig. 4. — Entropy per baryon as a function of height (z) from the center of the proto-neutron 
star. Solid line corresponds to that of Model 10CW. Sort-dashed line corresponds to that of 
Model 10C in QW96. 
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Fig. 5. — Outflow velocity, temperaure, density, and entropy per baryon as a function of 
height (z) from the center of the proto-neutron star in Model 10GS3. These quanta are 
written in unit of 10 7 cm s -1 , 1 MeV, 10 8 g cm -3 , and the boltzman constant, respectively. 
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Height (z) [km] 

Fig. 6.— Absolute values (in cgs units) of the components of Eqs. (16) and (17) for 
Model 10GS3. Lines (a)-(h) correspond to 4e/T, (P/ep)(e + P/p), v 2 z , pGM/z 2 , Pq/v z e, 
(e/ p + Pj p 2 )^, q/v z , and Q 2 D 2 /4 as a function of height (z) from the center of the proto- 
neutron star, respectively. The discontinuities of lines (e) and (g) at z ~ 8 x 10 6 cm and at 
z~2x 10 8 cm reflect the change of the sign of q and freezeout of the neutrino reactions 
(see subsection 2.1), respectively. 
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Table 1: Model names and input parameters. Mass and radius of the neutron star, total 
luminosity of neutrinos, mass outflow rate, and QD are showm respectively. 
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Table 2: Model names and output parameters. Entropy per baryon, dynamical 
timescase (r dyn ), analytically estimated dynamical timescale (Td yn . a na), electron fraction, final 
temperature are shown respectively. The reason why dynamical timescales are not written 
in some models is that temperature does not decrease to 0.2 MeV within z = 10 9 cm. 



